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Abstract 

We report on quantitative comparisons between simulation results of a bead-spring model and 
mode-coupling theory calculations for the structural and conformational dynamics of a supercooled, 
unentangled polymer melt. We find semiquantitative agreement between simulation and theory, 
except for processes that occur on intermediate length scales between the compressibility plateau 
and the amorphous halo of the static structure factor. Our results suggest that the onset of slow 
relaxation in a glass-forming melt can be described in terms of monomer-caging supplemented by 
chain connectivity. Furthermore, a unified atomistic description of glassy arrest and of conforma- 
tional fluctuations that (asymptotically) follow the Rouse model, emerges from our theory. 

PACS numbers: 61.25.Hq, 64.70.Pf, 61.20.Lc 
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I. INTRODUCTION 



Polymeric melts can often be cooled down easily to vitrify into disordered solids. It is an 
important challenge for first-principles approaches to develop an understanding of this tech- 
nologically important process. In polymer science a further important challenge is to derive 
well-known models of chain transport and relaxation in melts [l]. In a melt, excluded vol- 
ume interactions and chain connectivity cause subdiffusive segment motion, (approximately) 
described by the 'Rouse' and 'reptation' models which consider a single chain in an effective 
field pl. The challenge consists in deriving these models from microscopic interactions. 



In this paper, we propose an atomistic interpretation of the structura 



and conformational 
4j by quantitatively 



dynamics of a bead-spring model for an unentangled polymer melt 
comparing simulation and first-principles calculations. On the one hand, this explains the 
onset of the viscous slowing-down, ultimately leading to kinetic arrest into an amorphous 
solid (the glass transition) 5|, |6|, |7J, |8| . On the other hand, we find Rouse-like motion for very 
large chain length N, and also explain characteristic deviations from pure Rouse behavior 
for finite N. Our theory does not describe entanglements [^j because we start from isotropic, 
correlated monomer collisions which give rise to the "cage effect" in dense fluids, but vanish 
in the limit of infinitely thin chains where only topological constraints (entanglements) are 
present. We aim to describe fundamental consequences of the local steric packing in dense 
melts of flexible polymers and thus, in a first step, neglect chemical structure like torsional 
degrees of freedom; for simulation studies of the glass transition using chemically realistic 

nrin 

models see, e.g., Refs. [6l. 1 71 llOl]. 

Our approach is based on an extension of the mode-coupling theory (MCT) for the glass 



transition 



ll| to polymer systems. MCT predicts structural arrest - also referred to as the 



idealized liquid-glass transition - driven by the mutual blocking of a particle and its neighbors 
at a critical temperature T c which is located above the glass-transition temperature T g . 
Although complete structural arrest at T c is not observed in experiments and simulations, 

s properly 
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extensive tests of the theory carried out so far above T c suggest that MCT dea 
with some essential features of the structural relaxation in glass-forming liquids 

Our extension of the MCT to polymer systems will be done following the site formal- 
ism |l4|. In this formalism, each polymer molecule is divided into interaction sites, cor- 
responding to monomers or segments, and the dynamics as well as the static structure 
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of polymers are characterized by site-site correlation functions. For chains consisting of N 
monomers, the site-site correlation functions comprise 0(N 2 ) elements, and handling them is 
a formidable task for long chains. A key assumption of our 

_c theory Q 

the site-specific intermolecular surrounding of a monomer by an averaged one (equivalent- 
site approximation [ij], while keeping the full intramolecular site-dependence. For the 
statics this approximation has been verified by simulation for wave lengths around the av- 
erage segment separation for a bead-spring model of a polymer melt [l6j. Here we extend 
the test by presenting a quantitative comparison of the theory for collective and single-chain 
dynamics with molecular-dynamics (MD) simulations. This test requires as only input the 
average static structure factor S(q), the average intrachain structure factor w(q), and the 
site- resolved intrachain structure factors w^jq) to be defined below. These quantities are 



directly determined from the simulation 



161 ] . Such a fully quantitative comparison concern- 



ing the polymer dynamics is novel and, even for systems of simper constituents, has been 
done only recently [17]. It is motivated by successful MCT fits of the simulation data for 
our model (J). 

It might be appropriate at the start to summarize the philosophy of our approach for 
developing an atomistic theory of the slow structural dynamics in polymer melts, combining 
simulation and MCT. We aim to derive first-principles results for measurable quantities like 
intermediate scattering functions and mean-square displacements. Within MCT, such aver- 
aged information on the structural relaxation can be obtained from equilibrium structural 
input up to a single unknown parameter, the time scale parameter, which needs to be found 
from matching theory and simulation at one time instant (we will choose the final relax- 
ation time of density fluctuations at the wave length corresponding to the average segment 
separation for this purpose). We argue that quantities obtained by averaging, assuming 
homogeneity, should be studied first to characterize macromolecular motion, postponing 
consideration of heterogeneities and molecular/segmental variations to future more detailed 
studies. A central question of our investigation - besides the one whether such an approach 
is feasible at all - is which static information is required or sufficient to explain key features 
of macromolecular motion. Therefore, we consider a simple polymer model without chem- 
ical detail for which the required static input can be obtained with high precision and the 
theoretical predictions for the dynamics, calculated without adjustable parameters, can be 
compared with simulation results [17]. We expect to uncover fundamental mechanisms also 
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present in real polymer melts, for which, however, more complex static information than is 
necessary here might be required as input to a theory for the dynamics. 

The paper is organized as follows. Section [TT] introduces the model polymer system to 
be considered in the present article. In Sec. IIII[ we briefly review static properties of the 
model. Implications from the analysis of the static properties will be employed in Sec. IIVI 
to develop a tractable theory for polymer dynamics. In Sec. |V] theoretical predictions on 
structural and conformational dynamics of the model are compared with simulation results. 
The paper is summarized in Sec. |Vl]with some concluding remarks. Appendix |A] is devoted 
to the derivation of basic equations of motion, and Appendix [B] to the derivation of the 
Rouse model based on our microscopic approach. 



II. MODEL 



We studv a bead-spring model of linear chains, each containing N = 10 monomers of 



mass m [jj, U, ll8] . This is a model for highly flexible polymers, and is among the simplest 
models exhibiting glassy arrest and polymer specific dynamic anomalies. All monomers 
interact via a truncated and shifted Lennard- Jones (LJ) potential 

/4e LJ [(a LJ /r) 12 -(a LJ /r) 6 ]+C r<2r min , 
^Lj(r) = < (1) 
[0 r > 2 r min . 

In the following, all the quantities are expressed in LJ units with the unit of length olj, the 
unit of energy €lj (setting Boltzmann's constant ks = 1), and the unit of time (malj/e^j) 1 / 2 . 
The constant C = 127/4096 is chosen so that U^ir) vanishes continuously at r = 2r min 
with r min = 2 1 / 6 being the minimum position of the nontruncated potential. In addition, 
successive monomers in a chain interact via a finitely extendible nonlinear elastic (FENE) 
potential [3] 

h f r \ 

(2) 



1-1 ^ 



k 

^feneO) = --R 2 \n 

with Rq = 1.5 and k = 30. The superposition of the LJ and FENE potentials leads to a 
steep effective bond potential with a sharp minimum at r^, = 0.9606. 

For this model we carried out MD simulations of polymer melts at constant pressure p 
and constant temperature T. (The polymer melts comprise 100-120 chains, depending on 
temperature.) The MD simulations were performed in two steps js), 0]: For each T, the 
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volume of the simulation box was first determined in an isobaric simulation at p = 1. Then, 
this volume is kept fixed and the simulations are continued in the canonical ensemble using 
the Nose- Hoover thermostat. (The choice of this thermostat does not influence the relaxation 
dynamics of the melt jj], 8].) The simulations were carried out for the temperature range 
0.47 < T < 1, corresponding to monomer densities 0.91 < p m < 1.04. The lowest simulated 
temperature is slightly above T C MD m 0.45, the MCT critical temperature as determined from 
the MD simulation (cf. Sec. HVDj) . Before all measurements each state point (T, p m ) is fully 
equilibrated (the chains are allowed to diffuse several times over the distance corresponding 
to their radius of gyration). For each state point, all quantities are averaged over 150- 
200 independent time origins. A more detailed description of the simulation technique and 
simulation results for the model can be found in Refs. 3), 
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III. SUMMARY OF STATIC PROPERTIES 



Static structural and conformational properties of our model have been analyzed in detail 
in Ref. jl6|. In this section, we briefly review some of the main results of Ref. [if]], which 
help us to develop a tractable theory for polymer dynamics. 



A. Static structure factors 



Let us consider a polymer melt of n chains, consisting of N identical monomers, in 
a volume V. We denote by p = n/V the chain density. The static collective density 
fluctuations at the monomer level can most naturally be characterized by the monomer- 
monomer (or site-site) static structure factors 

S ab (q) = -<p a (q)*p 6 (q)>, (3) 
n 

defined in terms of the coherent monomer-density fluctuations for wave vector q 

n 

PM) =$>xp[2q-r 4 a ] (a = l,---,N). (4) 

i=i 

Here (•) denotes the canonical averaging for temperature T, and r° represents the position 
of the ath monomer in the ith chain. Since the melt is spatially homogeneous and isotropic, 
the structure factors depend only on the modulus of the wave vector, q = |q|. One can split 
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S a b(q) into intrachain and interchain parts 

S a b(q) = Wab{q) + phab(q), 



(5) 



in which the intrachain contribution is given by 

w <*® = ^E^Hq •(*?-*{)]), (6) 

and the interchain contribution by 

^(?) = -(E ex PHq-W-i?)] )• «") 

71 \ i+3 



These contributions reveal static correlations between monomers belonging to the same chain 
or to different chains, respectively. 

Commonly, not these site-resolved quantities, but structure factors averaged over all 
monomer pairs (a, b) are discussed. For instance, we obtain the collective structure factor 
of the melt by 

1 - 1 

S M = n E S M = ^(Ptot(q)Vtot(q)>, (8) 

a,b=l 

which can be represented in terms of the total monomer density fluctuations (the second 
equality in the above equation) 

N n N 

ptot(q) = E^(q) = EE exp [' q ' r ^]- ( 9 ) 

a=l i=l a=l 

The average S(q) can also be decomposed into intrachain and interchain parts 

S(q)=w(q)+p m h(q), (10) 

where p m = Np denotes the monomer density, and 

1 - 1 N 

w ( q } = iv E W M, %) = E h M- ( n ) 

a, 6=1 a, 6=1 



The intrachain contribution w(q) is often called 'form factor' in the polymer literature [2LI20]. 

Figure [1] shows the simulation result for S(q) of our model at T = 0.47, 0.70, and 1, which 
are representative temperatures in the investigated range 0.47 < T < 1.0. In this T interval, 
the collective structure of the melt is typical of a dense disordered system. Due to the weak 



compressibility of the melt, S(q) is small in the q — > limit. As q increases, S(q) increases 
toward a maximum which occurs around q* = 6.9 in our model. The corresponding length 
scale 2ir/q* is associated with the average segment separation which is of the order of the 
effective monomer diameter (= 1 within the reduced units). Thus, the dominant contribution 
to S(q*) comes from the amorphous packing in neighbor shells around a monomer. Upon 
lowering T, Fig. [1] indicates that the packing becomes tighter, which is reflected by the 
increased height of the peak S(q*) and by the shift of its position q* to larger values. Such 
changes of S(q) at q ~ q* reflect the interchain correlations since the intrachain structure 



factor w(q) was found to be nearly T independent [161 ]. 

In addition to monomer density fluctuations, the static spatial arrangement of the center 
of mass (CM) of chains and its T dependence might be of interest. The CM-CM static 
structure factor is defined by 



Sc(q) = -( J2 ex P["^ • ( R * " R ^ /' (12) 



n , 
uj=i 

where Rj denotes the CM position of the zth chain. The inset of Fig. [1] shows the simulation 
result for Sc(q)- It is seen that Sc(q) is fairly featureless: outside the small-g regime re- 
flecting the low compressibility of the melt, Sc(q) quickly approaches the ideal gas behavior, 
Sc(q) = 1- There is a tiny peak at qc ~ 3.4, whose height is, in contrast to that of S(q), 
practically T independent. 

One understands from Fig. [1] that the most pronounced T dependence in the static struc- 
ture occurs in S(q) around q* which reflect interchain monomer correlations. This implies 
that the slowing-down of the dynamics of our model upon lowering T cannot result from 
static CM-CM or intrachain correlations, but should be driven by interchain correlations 
at the monomer level, i.e., by the nearest neighbors that are not directly bonded to each 
other. As we will see, this is one of the principal predictions of our theory, according to 
which the coherent dynamics close to q* enslaves all other dynamics, including the CM and 
single-chain conformational dynamics. 

B. Equivalent-site approximation 

In Ref. [l^ . particular attention has been paid to the dependence of static correlation 
functions on the position of the monomer along the chain backbone to understand to what 



extent specific monomer-monomer correlations deviate from the average behavior. Since the 
site-site static correlation functions are necessary input quantities for the mode- coupling 
approach based on the site formalism, the comparison of the monomer- monomer correlation 
functions with their monomer-averaged counterparts can suggest suitable approximations 
and thus help developing a tractable theory. In this and the next subsections, we summarize 
such approximations which will be employed in Sec. HVl 

Let us introduce the site-site direct correlation function c a b(q) via the site-site Ornstein- 



Zernike equation 21| 



N 

Kb{q) = ^2 Cx y(i) + P h vb(q)]- ( 13 ) 

x,y=\ 

This is a generalized Ornstein-Zernike equation in which intrachain correlations are ac- 
counted for through w a b{q)\ it also serves as the defining equation of the direct correlation 
function in terms of w a b(q) and S a b(q), 

PCa b {q)=w~ b \q)-S- b \q). (14) 

Here X~ b (q) (X = w or S) denotes the (a, b) element of the inverse of the matrix X(g). 

The difficulty in dealing with the site-site correlation functions arises from the depen- 
dence on the indices (a, b). Such functions consisting of 0(N 2 ) elements cannot easily be 
handled for large N. One can argue, however, that, for long polymers, chain end effects for 
interchain correlation functions should be small, suggesting to treat all sites of a homopoly- 
mer equivalently. (This simplification is exact for a ring homopolymer.) This equivalent-site 
approximation is usually invoked for c a b(q), i-e., 

c(q) = c a b(q) (equivalent-site approximation). (15) 

Equation f|T5|) represents the principal idea of the PRISM (polymer reference interaction site 
model) theory developed by Schweizer, Curro, and coworkers [jjj]. 

Substituting Eq. ( fT5l) into Eq. ( fT3l) and then taking the summation b of the resulting 
equation, one gets the following scalar equation, called the PRISM equation, in terms of the 
averaged quantities defined in Eq. (fTT|) : 

h(q) = w(q)c(q)[w(q) + p m h(q)}. (16) 



S 



Equation (|T6|) provides the following expression for c(q) 



p m c(q) = l/w(q) - l/S(q), (17) 

in terms of the average w(q) and S(q). 

The validity of the equivalent-site approximation ffT5l) has been examined for our model 
by comparing c a b(q) obtained from Eq. (|T^|) with c(q) from Eq. (|T71) . with the quantities 
on the right-hand sides of these equations directly determined from simulations (see Fig. 5 



of Ref. 



16l|). It has been demonstrated that the approximation is well satisfied, except for 



functions involving the chain ends. This result suggests that, for our model, a theory for the 



melt dynamics can be derived by assuming Eq. ([15]) without introducing a large error 



22] 



C. Additional ring approximation 

Besides c a b(q), the static structure factors S a b(q) are necessary input quantities for the 
MCT based on the site formalism (cf. Sec. lIVj) . Thus, the equivalent-site approximation ffT5j) 
alone is insufficient to obtain a tractable theory since the specific monomer-position depen- 
dence still remains in S a b(q). This is obvious in view of the following relation 

S ab (q) = [{I - P w(g)c(g)}- 1 w(g)] ab , (18) 

which can be derived from Eqs. ([5]) and ( fi~3l) . Here I denotes the unit matrix. Thus, even with 
the assumption c a b(q) = c(g), a site dependence of S a b(q) results from chain connectivity, 
i.e., from the matrix structure of w a b(q). Therefore, it is desirable to have an additional 
approximation which simplifies the treatment of S a b(q). 

Remembering that the equivalent-site approximation (fT5j) is exact for a ring polymer, we 
will derive an additional approximation for linear chains based on another exact relation for 
rings. We will then examine the validity of this approximation in our simulation. 

A prominent feature of the site-site structure factor for a ring polymer is that S a (q) = 
J2b=iSab(<l) is independent of a, and there holds S a (q) = (1/N) X^Li ^a(o) — S(q); the 
second equality follows from Eq. ([8]). Furthermore, from the identity J2 X b S~^(q)S x b(q) = 1, 
we also have 5'~ 1 (g) = J2b=i ^ab (?) = V^*?)- m view of these exact relations for a ring 
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polymer, let us introduce the following approximations for linear chains: 

N 

£ a (g) = X^<?) ~ (19) 

6=1 

N 1 

S^^J^^^Siq)- (20) 

Figure [2] and its inset examine to what extent the ring approximations (fT9l) and fl20l) hold 
for our model. It is seen that, except for Si(q) and S'f 1 (g) referring to the end monomer, 
the ring approximation is well satisfied, suggesting that this additional approximation can 
also be used in deriving a theory for polymer dynamics without introducing a large error. 
Let us add that both the equivalent-site and ring approximations have been found to hold 



well also for the liquid structure of some semiflexible polymer models 



23|. 



IV. THEORY 



A. MCT equations for coherent structural dynamics 

In the site formalism, collective structural dynamics are to be described by site-site density 
correlators 

F ab (q,t) = -(p a (qYe lCt p b (q)) (a, b = 1, ■ ■ ■ , N), (21) 
n 

whose initial values are the static structure factors S a b(q) = F ab (q,0). Here, C denotes 
the Liouville operator appropriate for Newtonian dynamics. MCT equations of motion 
for F ab (q,t) for general flexible molecules are derived in Appendix |Aj and consist of the 
Zwanzig-Mori exact equation of motion and an approximate expression for the memory 
kernel. The former is obtained by introducing a projection operator V onto the monomer- 
density fluctuations and the corresponding longitudinal current fluctuations, and reads [cf. 
Eq. (|A6|)]: 

N N t 

d*F ab {q,t)+J2nlM F xb{q,t) + J2 / dt'M ax (q,t-t')d t ,F xb {q,t') = 0. (22) 
x=l «=i Jo 

Here VL 2 ab (q) represents the characteristic frequency given by 

n^(ff) = gVS#( g ), (23) 
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with v 2 = k^T jm (= T within the reduced units) denoting the monomer thermal velocity, 
and a formally exact expression for the memory kernel reads 

M ab (q,t) = -^(f a (q)*ex V {iQCQt)f b (q)), (24) 

in terms of the fluctuating random force / a (q) which evolves with the generator QCQ, 
where Q = 1 — V. An approximate expression for / a (q) as derived under the mode-coupling 
approach in Appendix |A] is given by [omitting the irrelevant factor — % from Eq. ( 1A211) 

2 N 

f*w = vE E(q • k ) p*o*) ^(p)' ( 25 ) 

k x=l 

in which q = q/g and p = q — k. With the use of the factorization approximation (]A13I) . 
one finally arrives at the following MCT expression for the kernel [cf. Eq. (1A22I) ]: 

2 N r 

M ab (q,t) = / dk{(q-k) 2 c ax (k)c by (k)F xy (k,t)F ab (p,t) + 

^ ' x,y=l 

+ (q-k)(q-p)c ax (k)c by (p)F xb (k,t)F ay (p,t) } . (26) 

Equations (1221) and (1261) provide a set of closed equations for determining site-site coherent 
density correlators F ab (q,t), provided static quantities S ab (q) and c ab (q) are known. From a 
computational point of view, however, it is quite demanding to solve these N x N matrix 



equations since iV may become large for polymeric systems [241 ] . It is at this point where 
the analysis of the static properties, presented in Ref. and summarized in Sec. IHIl will 
help us to develop further approximations. 

As mentioned in Sec. IIII Bl the equivalent-site approximation, c ab (q) = c(q), is well 
justified for our model. So, we insert c ab (q) = c(q) in Eq. (1251) and obtain 

/.(q) = — £(q ■ k) c(k) Ptot (k) Pa ( P ), (27) 



n 

k 



where ptot(k) denotes the total monomer density fluctuations introduced in Eq. (Q. This 
expression reveals that the equivalent-site approximation alone does not suffice to simplify 
the problem: the dependence of / a (q) on the monomer position remains, and the resulting 
M ab (q,t) still carries the (a, h) dependence, i.e., it consists of 0(N 2 ) elements. Furthermore, 
no simplification is yet achieved concerning the frequency matrix (l2"5j) . 

Progress is made if we invoke the second approximation described in Sec. IIII CI The 
frequency term (1231) can be simplified by the use of the ring approximation (120]) to 



N N 



E^(^) = ? V E^) « l^lSiq) = V 2 (q). (28) 



a=l a=l 
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A corresponding simplification can be introduced for the fluctuating force / a (q). To this 
end, we notice that / (q) in Eq. fl271) originates from interchain interactions represented 
through the direct correlation function. It is then reasonable to expect that, for long chains, 
the specific monomer-position dependence in / a (q) is small, and can be well approximated 
by an averaged one. (This approximation is exact for rings.) We therefore introduce an 
approximation 

1 N 

/a(q)«^£/a(q)- (29) 

a=l 

This ring approximation replaces the site-specific surroundings of a monomer by an averaged 
one, and then the fluctuating force is given by 



■kfa) = ^ E@ ' k ) C ^ Ptot(p). (30) 

k 

This leads to an expression for M a b(q,t) (to be summarized below) which now does not 
depend on the site indices (a, b). 

The approximations discussed so far allow us to derive a set of closed MCT equations for 
the collective total monomer density correlators 

1 - 1 

F fa*) = N S = ^(Aot(q)*e' £ Vtot(q)), (31) 

o,6=l 

whose initial value is F(q, 0) = S{q) [cf. Eq. ©]. To this end, we take (1/JV) J2 a ,b of E 1- fl2D> 
and then insert the frequency term (|28|) and the memory kernel with the fluctuating force 
given in Eq. (I3"U1) under the factorization approximation (IA13I) . This gives the following set 



of MCT equations for the normalized coherent density correlators 4>(q,t) = F(q,t)/ S(q): 

d^(q,t) + Q 2 (q)<P(q,t) + Q 2 (q) [ dt' m(q,t - if) d t ,<j>(q,t') = 0, (32) 

J o 

rn{q, t) = ^Jdk V(q; k, p) <f>{k, t) <f>{p, t). (33) 

Here Q 2 {q) = q 2 v 2 /S(q), and the vertex function reads 

V(q; k, p) = S(q)S(k)S(p) {q • [kc(k) + pc(p)}} 2 . (34) 

One can solve these equations for cj)(q, t) provided the average static quantities S(q) and c(q) 
are given as input. 

Equations (1321) and (|33|) merit some comments: (i) These equations are formally identi- 
cal to MCT equations for monatomic liquids. Polymer-specific effects, such as local stiffness 

12 



of the chain backbone or chain length N, enter the relaxation only via the direct correla- 
tion function c(q), the structure factor S(q), and the monomer density p m . These static 
equilibrium features fully determine the long-time coherent dynamics of the melt, (ii) Equa- 
tions f[3"31 and (|34[) indicate that the memory kernel contains the factor S (q) S (k) S (p) . So, 
the slow dynamics upon lowering T should be mainly driven by wave vectors close to q* be- 
cause there, S(q) is largest and the strongest dependence on T occurs (c/. Sec. IIII Bj) . Thus, 
our theory predicts that the glassy structural slowing-down is connected to the increase of 
the first peak of S(q), i.e., to the local cage effect. 



We finally notice that the so-called regular contribution to the memory kernel [U[ is 
discarded in our theory, and our approximate memory kernel is completely given by the 
mode-coupling expression. The latter provides the slow contribution relevant for the struc- 
tural slowing-down. The regular contribution is supposed to embody memory effects already 
present in the normal high-T state of liquids, and accounts for the fast dynamics in the short- 
time regime. We drop the regular contribution since it does not affect the MCT predictions 
on the slow- relaxation regime III, [25] . Thus, care has to be taken in comparing theoretical 
predictions with simulation results, since the theory without the regular contribution does 



not properly describe the short-time dynamics 



26|. 



B. MCT equations for single-chain dynamics 

The basic variable characterizing the dynamics of a single (or tagged) chain is 

Pli^t) = e^W, (35) 

where r"(i) denotes the position of ath monomer in the tagged (labeled s) chain at time t. 
The density correlator for the single-chain dynamics is defined by 

F: b (q,t) = (pl( q ypl(q,t)), (36) 

whose initial value is the intrachain structure factor, w a b(q) = F^ b (q, 0). 

The derivation of the MCT equations for F* b (q,t) is outlined in Appendix IA 31 and the 
resulting matrix equations can be summarized as 

N N t 

d?F: b (q,t) + Y,nil(q)F: b (q,t)+J2 nS ^ I dt f ml y (q,t - t') d t ,F* b (q,t') = 0, (37) 

x=l x,y=l 

13 



where the frequency matrix is given by 

n s a Uq) = q 2 v 2 w- b 1 (q) 1 (38) 

and the MCT expression for m s ab (q,t) under the equivalent-site approximation (1151) reads 

- f 

m s ab (q,t) = J2^Aq) / dk V s (q; k, p) F^k, t) <j>{p, t), (39) 

x=l J 

with the vertex function 

Hq; k, P ) = j^s{ P ) (q ■ P ) 2 c{ P f. (40) 

Equations (}3~Tj) and (|39|) constitute a set of closed N x iV-matrix MCT equations for the 
single-chain density correlators F* b (q,t). One can solve these equations with the knowledge 
of the static quantities - S(q), c(q), and w a b(q) - and of the coherent density correlators 
4>(q, t). It is clear from Eq. ( |39l) that the slowing-down of the single-chain dynamics is driven 
by that of the coherent dynamics. 

Unlike the MCT equations for the coherent dynamics, one cannot simplify the matrix 
structure of Eqs. (13T1) - (13"9"1) for the single-chain density correlators F^ b (q, t) in order to prop- 
erly describe chain-connectivity effects, taken into account through the intrachain structure 
factor matrix w a b(q). For example, one needs the site-site F^ b (q, t) to fully describe the chain 
conformational dynamics, i.e., all the Rouse- mode correlators introduced below, and has to 
solve the matrix MCT equations ( 1371) and ( |39i) for this purpose. Let us note in this connec- 
tion that, from a computational point of view, it is not so demanding to solve these matrix 
MCT equations for F* b (q,t). This is because the most time consuming part in numerically 



solving the MCT equations is spent in solving the ones for coherent dynamics 



27] 



C. MCT equations for Rouse-mode correlators 

In this subsection, the site-density description of the single-chain dynamics shall be related 
to the traditional Rouse description, and we derive the MCT equations for the Rouse-mode 
correlators. Let us stress that the present rewriting is exact and always possible. The N 
degrees of freedom of segmental motion are mapped onto N modes labeled by p. If the 
Rouse model holds, the modes will be statistically independent and the matrix of Rouse- 
mode correlators introduced below becomes diagonal. 
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Let us introduce the following N x N orthogonal matrix P ap (a = 1, 2, • • • , N and p 
0,1,..- ,N-1): 



ap 



2 



cos 



(a - l/2)pn 
N 



(p = 0), 

(p=l,2,.-. ,AT-1) 



(41) 



The Rouse-mode vectors shall be defined in terms of the monomer positions as X p (t) 

Za=lPap<(t) |2a|,i.e., 



r rr * 

o=l 

A 7 " 



cos 



(a - l/2)jwr 



(p = 0), 



(p = i,2,... ,iv-i; 



(42) 



The inverse relation is given by r" (t) = J2 p= q -Pa P X p (t): 



JV-l 



r:(t) = V ^ 1 Xo(t) + ^ E x ^) cos 



p=l 



(a - l/2)p?r 
N 



(43) 



We introduce the Rouse-mode correlators as 



[(X (0) • X (i)> - (X (0) • X (0))] / 3N (p = p' = 0), 

C pp ,{t) = { [(X (0) ■ Xp,(t)> - (X (0) • X p ,(0)>] JZN (p = 0,p'^ 0), 

(X p (0)-X p ,(t))/3iV (P^O.^O). 

From the definition, it is obvious that 

Coo(0) = 0, C 0p (0) = for p > 0. 



(44) 



(45) 



Let us introduce the following (N — 1) x (iV — 1) matrix to denote the initial values of C pp '(t) 
for p,p' > 0: 

C pp i = Cpp'(0) defined only for p,p' > 0. (46) 

It is necessary to introduce this new matrix in order to discuss the inverse matrix of C pp >(0): 
the inverse of the N x N matrix C pp i(0) does not exist because of Eq. (j45D . whereas that of 
the (N — 1) x (N — 1) matrix C pp / does exist, and shall be denoted as C~2. 

Since X (i) = y/NH s (t) with R s (i) denoting the CM position of the tagged chain at time 
t, Coo(i) is related to the CM mean-square displacement (MSD) gc(t) = ([R s (t) — R s (0)] 2 ): 



g c {t) = -6C 00 (t). 



(47) 
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For the monomer MSD, g a (t) = ([r°(t) - if (0)] 2 ), it follows from Eqs. (035 and (JED that 

■ V_J r (a-l/2)p'n 



N 



g a (t) = gc(t) - 12V2 Y^ c wW) cos 

p'=l 

N-l JV-1 

+ PpA®) - Crf(t)} cos 

p=l p'=l 



(a - l/2)jwr 
iV 



cos 



(a - l/2)j/7T 
N 



(48) 



In particular, one gets for the MSD averaged over all the monomers in a chain 

N N-l 

g M (t) = (l/N)J2da(t) = 9c(t) + 6 ^[C pp (0) - C„(t) }. 

a=l p=l 



(49) 



We next consider how the Rouse- mode correlators C pp i{t) are related to the site-density 
correlators F* b (q,t). This will allow us to write down MCT equations for the former based 
on the ones for the latter. Since the density fluctuations of the tagged chain for small q are 
given by p s a (q,t) ~ 1 + zq • T*(t), it is easily understood that C pp >(t) can be expressed as a 
linear combination of F* b (q,t) for q — > 0. Indeed, one can show that 



1 N 



o,b=l 



l+q 2 [C 00 (t)-A}+O(q 4 ) (p = P ' = 0), 

q 2 [C Qp ,(t) ~ B pl ] + 0(g 4 ) (p = 0,p'> 0), (50) 



q 2 C pp ,{t)+0{q 4 



ip,p' > 0), 



where A = (1/6 A/" ) X]a6=i(( r s ~~ r s) ) = -Rg/3 with i? g denoting the radius of gyration of a 
chain, and B p = (1/3\^2N 2 ) J2ab=i(( r t~ r l) 2 ) cos[(b — l/2)pTt/N). Using these relations, one 
derives the following MCT equations for C pp /(t) by taking the q — > limit of Eqs. (137 1) -( 1401) : 

N-l N-l „ t 

d 2 C pp ,(t) + v 2 D pp , + v 2 E pp „C p> ,(t) +v 2 Y dt'm pp// (t - t') d t ,C p n p ,{t') = 0, (51) 

^// — n ^11 — n 



p"=0' 



where 



Dppi — 5 p0 5 p i /N, 



E pp > 



c; p )/n (p,j/>o), 

(otherwise). 



(52) 
(53) 



and the expression for the memory kernel reads 



m pp/ (t) 



67T 2 



dk k 4 S(k)c(kf 



N 



PapF^k,^ 



bp' 



.a, 6=1 



<f>(k,t). 



(54) 
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These MCT equations for C pp r(t) can be solved provided the static quantities - S(q), c(q), 
w a b(q), and C vv > - and the full site-site single-chain density correlators F^ b (q, t) as well as the 
coherent density correlators <j)(q, t) are known. MSDs can then be obtained from Eqs. (I47p 
through (|49p . Again, it is clear from Eq. (|54|) that the slowing-down of the Rouse-mode 
dynamics and MSDs is dictated by that of the coherent dynamics. 

In Appendix [Bj we will show how the Rouse model emerges from our MCT equations in 
the asymptotic limit of large N. 



D. Universal MCT predictions 

Here, we briefly summarize some universal MCT predictions which are necessary for 
understanding the present paper. As described in Ref. 29], all universal results concerning 
the MCT-liquid-glass transition dynamics, originally developed for simple systems 111 ], are 
also valid for molecular systems, and the MCT for polymer melts developed in Sec. II VI shares 
this feature. This justifies the use of the MCT universal predictions in analyzing polymer 



data, whose validity has been tacitly assumed in previous studies of our model \4, la. |3C ] . 

One of the central predictions of MCT is the existence of a critical temperature T c . 
The long-time limit, or the nonergodicity parameter, of the coherent density correlator 
f(q) = 4>(q,t — > oo) obeys the implicit equation 

d%)=*M. («) 

which can be derived by taking the t — > oo limit of Eqs. ( 1321) and ( 1331) and introducing the 
mode-coupling functional Fqlf] = (1/2) / rfkV(q; k, p)f(k)f(p). One gets trivial solutions 
f(q) =0 for T > T c , meaning that the density fluctuations relax completely at long times, 
a characteristic feature of the ergodic liquid state. On the other hand, nontrivial solutions 
f(q) > can be obtained for T < T c , describing nonergodic dynamics in which density 
fluctuations cannot fully decay. The nonergodicity parameter f(q) measures the "solidity" 
of such an amorphous solid on length scales ~ 2n/q, and is thus also referred to as the 
glass-form factor or the Debye- Waller factor. The ergodic-to-nonergodic transition at T c is 
called the idealized glass transition, and f(q) at T = T c , to be denoted as f°(q), is referred 
to as the critical nonergodicity parameter. It also has the meaning of the plateau height in 
the two-step relaxation of <p(q, t), and quantifies the strength of its a relaxation (see below). 
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MCT predicts that for temperatures close to but above T c , to which we restrict our 
attention in the following, the dynamics of any time- dependent correlation function <px(f) 
coupling to density fluctuations exhibits a two-step relaxation: the relaxation toward the 
plateau, and the final relaxation from the plateau to zero. These relaxations are respectively 
characterized by the time scales t a and t' a defined by 

t a = t /\a\ s , 5 = l/2a, (56) 
t' a = t B-^ b 1 1 a\\ 7 = (l/2a) + (1/26). (57) 



Here to denotes some microscopic time scale [311 ] . a is called the separation parameter, and 
measures the distance from the critical point: 

a = Ce, e=(T c -T)/T c . (58) 

Except for to, all the exponents and the constants B and C appearing in these equations can 



be evaluated from the mode-coupling functional J~ q [f] with the knowledge of f c (q) [ll|, 132 ]. 
According to MCT, the dynamics dramatically slows down upon lowering T since the time 
scales t a and t' a diverge for T — > T c +. The two-step-relaxation scenario emerges since the T 
dependence of t' a is stronger than that of t a . 

The dynamics which occurs near the plateau is referred to as the j3 process. The height 
of the plateau is given by the critical nonergodicity parameter f x of the correlator <px{t), 
which can be determined from the corresponding mode-coupling functional via an equation 
analogous to Eq. (I5"5"j) . MCT predicts that there holds for o — > — (i.e., for T — > T c +) 

Mt) = fl + h xv ^\g(t/U) + 0(a). (59) 



Here hx is called the critical amplitude and g(i) the (3 correlator H[. Equation ( |59l) is called 
the factorization theorem, according to which the dependence of the correlator <px{t) on X 
(e.g., the wave number) represented through hx is factored from the temperature and time 
dependence described by ^/\a~\g(t/t a ). 

The decay of 4>x{t) down from the plateau f x is called the a process. For this process, 
MCT predicts for o — >• — 

Mt) = Mt/l), (60) 

which is also referred to as the superposition principle. The temperature independent shape 
function <px{t) - referred to as the a-master function - is to be evaluated from the MCT 
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equations at T = T c , and the temperature dependence is given via the time scale t' a intro- 
duced in Eq. fl57|) . The initial part of the a process is given by the von-Schweidler law, 
(fix(t) = fx — hx(t/t' a ) b [ill ]. The superposition principle implies that the a-relaxation time 
Tx of any dynamical variable, defined e.g. via the convention </>x( r x) = 0.1 which will be 
employed in this paper, is proportional to t' a , i.e., 

t x = CxC (61) 

This implies that the unspecified microscopic time scale to, which is common to all dynamical 
quantities, can be eliminated by choosing a particular variable X, e.g., the collective density 
correlator <f)(q,t) at the structure factor peak position q*, and then plotting other quantities 
as a function of t/r q * where r q * denotes the a-relaxation time of (p(q*,t). 

For our model, the mentioned MCT universal predictions have been successfully ap- 
plied to analyze simulation data, from which various characteristic quantities have been 



extracted 



4, 



30] . For example, the critical temperature T C MD « 0.45 was obtained from 



the consistent analysis of both (3 and a relaxations. Critical nonergodicity parameters have 
been determined for several correlators by applying Eq. (j59p including its leading correc- 

ri 

tion [32] in the (3 regime. These simulation results can directly be compared with our new 
first-principles theoretical calculations. 

In addition, and more importantly, our microscopic theory developed here can make 
predictions concerning polymer specific features which are outside the scope of the univer- 
sal MCT predictions. Through the comparative study of such first-principles theoretical 
predictions and simulation results, we will propose an atomistic interpretation of the slow 
structural and conformational dynamics of supercooled polymer melts. 

V. RESULTS AND DISCUSSIONS 
A. Collective structural dynamics 

We start by comparing theoretical and simulation results for the critical glass-form factors 
f c (q) of the coherent density correlators <j)(q,t). The theoretical result for f c (q) can be 
obtained by solving Eq. (1551) provided static inputs at T C MCT are known. (T c determined 



from the theory based on the analysis of Eq. 



551) will be denoted as T MCT to discriminate it 



from T C MD obtained from previous analyses 4, 8( of the simulation data.) If T C MCT lies in the 
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range of temperatures for which simulation results are available, the static inputs at T^ 1 ^ 
can be determined accurately. However, a problem arises if T C MCT is found to be below the 
lowest simulated T. This was the case in our study. (See Sec. I VI I for a possible explanation.) 
So, we had to estimate the required static input by a linear extrapolation based on the 
simulated S(q) at T = 0.47 and 0.48. (No extrapolation was necessary for static intrachain 
correlation functions since they are nearly independent of T as shown in Ref. [li]].) The 
dashed line in Fig [1(a) shows the extrapolated S{q) at T C MCT » 0.277 (< T C MD « 0.45) and 
the inset compares its behavior close to the first-peak position q* with the simulated S(q) 
at T = 0.47, 0.48 and 1. For the simulated S(q) the peak height increases and the peak 
position shifts to larger q upon lowering T as discussed in Sec. IHI Al The extrapolated S(q) 
inherits this trend. This suggests that the physics should not be significantly altered due to 
possible errors in our extrapolation. 

Figure EJ^a) compares f c (q) from MCT with the simulation result determined in Ref. {^]. 
For q > q* = 6.9, i.e., for distances comparable to the average monomer separation, we 
find a high degree of accord between theory and simulation. In particular, the agreement is 
quite good at q*, which is gratifying since the coherent dynamics for wave vectors close to q* 
drives the glassy slowing-down [cf. the second comment below Eq. (j34|) ] . On the other hand, 
the theory fails to reproduce the shoulder present in the simulation results at intermediate 
q near the peak qc = 3.4 of Sc{q)- 

A similar conclusion can be drawn for the a-relaxation time r q , defined via the convention 
4>(q,T q ) = 0.1, which is shown in Fig. [3]^b). r q depends nonmonotonically on q and varies 
over an order of magnitude from r q * « 1500 to r g= i6 ~ 40 [4]. The theory semiquantitatively 
captures these trends for q > q*, but misses the peak around qc- We also find corresponding 
deviations at q ~ qc in the relaxation stretching, as demonstrated in the inset of Fig. [3](b), 
quantified in terms of the stretching exponent (3 q which is obtained via a Kohlrausch-law fit 
of the a-decay part of the correlator, <p(q,t) oc exp[— {t/r')^]. 

Circles in Fig. H|a) shows 4>(q,t) for q = 4.0, 6.9, and 12.8, obtained from simulations 
at T = 0.47 which is close to T C MD rs 0.45. Clear evidence for the presence of a two- 
step relaxation exists. As described in Sec. IIVD| MCT provides, up to a time scale to 
common to all dynamical quantities, quantitative predictions for the a regime in terms of 
the (T independent) a-master curves, which are drawn as solid lines in Fig. Ufa). The 
dynamics including the early (3 regime can be described by solving the MCT equations 
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for a temperature T above T C MCT , which will be characterized by the distance parameter 
6 mct _ ^jtmct _ 7 1 ^t i mct_ rj^g dashed curves in Fig. 11(a) present such theoretical results 
for e MCT = —0.046 which corresponds to the distance between T = 0.47 and T C MD pa 0.45 of 
the simulation result. Since the separation parameter controls the ratio of the time scales 
characterizing the (3- and a-relaxation regimes [cf. Eqs. (1561) and (157)) ]. we found that a 
better agreement in the early f3 regime can be achieved by treating e MCT as a fit parameter. 
The theoretical results for e MCT = —0.022, which has been chosen so as to reproduce better 
the time-scale ratio found in the simulation data, are shown as dotted lines in Fig. ID(b). 
All the theoretical and simulation curves shown in Figs. |D(a) and IU(b) are plotted versus 
t/r q * with r q * being the a-relaxation time at q*. In this way, the unspecified time scale to 
can be eliminated [cf. discussion below Eq. (151))]. and the theoretical prediction on the q 
dependence of the a-relaxation times can be examined. 

Figure Hla) and 0](b) demonstrate that, for q = 6.9 and 12.8, the MCT curves quanti- 
tatively describe the simulation results in both the /3 and a regimes including the relative 
magnitude of r q and the stretching of the relaxation. (The disagreement for short times 
arises mainly because, as mentioned at the end of Sec. IIV Al the regular part of the memory 
kernel is not included in our theory.) On the other hand, the agreement is not satisfactory 
at q pa q c . We thus conclude from Figs. [3] and H) that, except for a process that occurs 
at q pa g c , our theory describes the coherent structural dynamics of polymer melts at a 
semiquantitative level. Since, as will be discussed at the end of Sec. IV Dl the dynamics at 
q pa qc does not appear to be directly related to the glass transition, this verifies one of the 
principal predictions of our theory: the emergence of the glassy slow dynamics is connected 
to the increase of the first peak of S(q), i.e., to the local cage effect. We will further comment 
on the dynamics at q pa q c below. 

B. Single-chain density correlators 

We next consider the single-chain density correlators. This will be done in terms of 
the averaged single segment correlators <f) s (q,t) and the (normalized) collective single chain 
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correlators <fi p (q,t) defined by 



N 



(P s (q,t) 



(W)^F a s a (g,t), 



(62) 



a=l 



N 



4?{q,t) 



(1/iV) Mff)- 



(63) 



a,b=l 



4> s (q,t) is denned solely in terms of the diagonal (a = b) elements of F^ b (q,t), i.e., it probes 
only the self-motion of monomers, whereas <ft p (q,t) reflects also interference effects from 
other monomers belonging to the same chain through the off-diagonal (a ^ b) elements. 
The simulation results for these correlators for our model have been analyzed in Ref. {4]. 

Theoretical and simulation results for the critical nonergodicity parameters of <j) s (q,t) 
and (f) p (q,t), to be denoted as f sc (q) and f pc (q), are compared in Fig. [5(a) and [5(b). The 
theoretical results are determined from the equation analogous to Eq. ( 1551) . which can be 
obtained by taking the t — > 00 limit of Eqs. (157)) and (1591) . The full time-dependence of 
the correlators <p s (q,t) and <p p (q,t) for representative wave numbers are shown in Figs. [6(a) 
and [6(b), where the MCT a-master curves and the MCT curves at the distance parameter 
e T = —0.046 are compared with the simulated correlators at T = 0.47. All the curves 
in Figs. [6(a) and [6(b) are plotted versus t/r q * with r q * being the a-relaxation time of the 
coherent density correlator <ft(q,t) at q = q*. In this way, the theoretical predictions not 
only on the q dependence of the a-relaxation times of 4> s (q,t) and 4> p (q,t), but also on the 
relative time scale of the single-chain and coherent dynamics can be tested. 

We first notice from Fig. [3(a) that the theory (solid line) describes well the simulated 
f sc (q) (circles) including the wave numbers q < q*. Apparently, this seems inconsistent 
with the result for the coherent dynamics where we found disagreement at q w qc, but 
can be understood in the following way. In the small-g limit, there holds <fi s (q — > 0,t) = 
1 — ( ( ? 2 /6)5'm(^) +0(q 4 ) in terms of the monomer averaged MSD, which implies the Gaussian 
approximation (p s G (q,t) ~ exp[— (q 2 /6)gM(t)}. The discussion of gu{t) will be presented in 
Sec. IV Dl but let us mention here that the critical nonergodicity parameter 6(r^) 2 of <?m(£) 
- the plateau height of #m(£) in the j3 regime - quantifies the size of the cage of a monomer 
formed by its surroundings. Thus, one expects the approximation 



to be valid at least for small q. The dash-dotted line in Fig. [5](a) shows fs c (q) with the value 



/-(g)^exp[-g 2 (r^) 2 ], 



(64) 
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= 0.098 taken from the theoretical calculation. It is seen that fQ°{q) describes well both 
the theoretical and simulated f sc (q). This reveals that the g-dependence of f sc (q) mostly 
reflects the spatial extent of the localized monomer motion in the /3-relaxation regime, and 
that our theory quantitatively describes the simulated f sc {q), since it correctly predicts the 
size of the cage. Figure 0(a) furthermore indicates that, not only the dynamics in the f3- 
relaxation regime where the correlators <fi s (q, t) are close to the plateaus f sc {q), but also their 
a decay is well described by the theory, including the relative time scale of the single-chain 
and coherent dynamics. 

In contrast to f sc {q), we expect some interference effects due to chain connectivity to be 
visible in f pc (q) on top of the nearly Gaussian background just mentioned. The theoretical 
prediction for f pc (q), shown as the solid line in Fig. G^b), indeed exhibits an oscillatory 
feature which is in phase with w(q). The oscillatory q variation of w(q) in turn reflects 
the bonding effect as discussed in Ref. In the simulated f pc (q) } shown as circles in 

Fig. E](b), the presence of such oscillations as predicted by the theory is discernible, though 
its amplitude is much smaller. (The oscillation in the simulated f pc {q), though tiny, can 
more easily be grasped in Fig. 8 of the first article in Ref. Q] where both of the simulated 
f sc (q) and f pc (q) are plotted in one panel.) This explains, e.g., why the theory does not 
describe so well the simulated (p p (q,t) at q — 4.0 [Fig. E](b)], in spite of the fact that it 
well describes the simulated <fr s (q,t) at the same wave number [Fig. E^a)]: the oscillatory 
q dependence of the plateau height and the a-relaxation time is more pronounced in the 
theoretical prediction. 

C. Rouse-mode dynamics 

We next turn our attention to the Rouse-mode correlators C pp >(t) describing the chain 
conformational dynamics. Circles in Fig. [3(a) show the simulation results for the normalized 
Rouse-mode correlators c p (t) = C pp (t) / C pp (0) at T = 0.47 for representative mode indices 
p. It is seen that c p (t) do not clearly exhibit the two-step relaxation. This is because the 
plateaus /° of c p (t) are so large [/£ > 0.9 for simulation results as shown in Fig. [11(a)] that 
only about or less than 10% of the decay is left for the relaxation towards the plateau. 
Thus, most of the relaxation of c p (t) occurs in the a regime. We therefore included for 
comparison only the MCT a-master curves as solid lines in Fig. 0(a), and the comparison 
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including the early (3 regime is done separately in Fig. [3(b) which highlights the plateau 
regime. For a more quantitative comparison of the features in the a relaxation, theoretical 
and simulation results for the plateau heights f£, the a-relaxation time r p defined via the 
convention c p (r p ) = 0.1, and the stretching exponent /3 P based on the Kohlrausch-law fit of 
c p (t) are compared in Fig. [S]for all the Rouse modes p. 

It is seen from Figs. [7] and [8] that our first-principles theory describes at a semiquantita- 
tive level main features of the simulation results for c p (t), such as the high plateau values f£ 
and the nearly exponential relaxation (/3 p > 0.9). In particular, from the comparison of the 
ratio T p /r q * of the a-relaxation times shown in Fig. [51(b) , we see that the theory provides 
a good description of the time-scale separation of the single-chain conformational fluctua- 
tions (characterized by r p ) from the local dynamics of the surrounding medium (r g »), which 
becomes more pronounced with decreasing p. 

The Rouse theory assumes a chain to be in a Markovian heat bath, i.e., that all dy- 
namical correlations in the surroundings are much faster than the single-chain dynamics [2[ . 
However, since a polymer is surrounded by identical polymers, the assumption of the time- 
scale separation cannot be justified a priori. Our microscopic theory developed here verifies 
this central assumption from first principles: our theory predicts the time-scale separation 
for small p/N, and hence, the Markovian approximation can be justified for small Rouse- 
mode indices (see Appendix [B]). The simulated results for the Rouse-mode correlators do 
not exhibit pure Rouse behavior due to finite- N effects (cf. Appendix IB 4[) . and even such 
deviations can be semiquantitatively accounted for by our theory as demonstrated in Fig. M 
where the pure Rouse behavior is also included for comparison. 

On the other hand, we also observe features predicted by the Rouse model, both in the 
simulated and theoretical results, already for such short chains as iV = 10. Let us recall 
from Ref. 33[] that the matrix of the equilibrium values of the (unnormalized) Rouse-mode 



correlators for our model is close to diagonal, C pp '(t = 0) ~ S pp f, and that the amplitude 
C pp (t = 0) for the smallest p is approximately given by the Gaussian result. It is an important 
numerical observation within our approach that the Rouse- mode correlators C pp >(t) remain 
nearly diagonal for all the times, in agreement with the simulation result, and that the 
time scale separation holds rather well between the memory functions in Eq. (I54j) and the 
Rouse-mode correlators for the smallest p. This holds because collective density fluctuations 
at microscopic wave lengths dominate the memory functions. This is also the reason why 
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the theory yields asymptotically in the large-iV limit the Rouse model spectrum and the 
characteristic y/i anomaly in the average segmental MSD (see Appendix IB1). 



D. Mean-square displacements 

Now, let us see how the Rouse-mode dynamics affects the single-chain diffusional pro- 
cesses. Circles in Fig. [9]shows the simulation results for the monomer-averaged MSD, gwi{t), 
and the CM MSD, gc(t). The MSDs also exhibit a two-step relaxation: after the short-time 
ballistic regime gxif) oc t 2 (X = M or C), the increase of the MSD begins to be suppressed 
due to the cage effect, and there appears the (3 regime where gxif) is close to a plateau 
which will be denoted as 6(r^) 2 . The appearance of the plateau regime reflects the confined 
dynamics of monomers inside the cage, and the height of the plateau of the monomer MSD 
gwi(t) reflects the size of the cage. We find both from theory and simulation « 0.1 (cf. 
Fig. [9]), i.e., the amplitude of the confined dynamics inside the cage is about 10% of the 
monomer diameter. The increase of gx(t) above the plateau towards the diffusion asymp- 
tote, gx(t) = QDt with D denoting the diffusion constant, is the a process of the MSD. In 
contrast to gc(t), gu(t) i n this regime is significantly affected by chain connectivity since 
the monomers participate in the conformational motion and most of the conformational 
fluctuations reflected by c p (t) occurs in the a regime as mentioned above. As a result, a 
polymer specific anomaly in the MSD - a subdiffusive (~ t x ) regime - emerges in gyi(t) 
in the a regime. Since MCT predicts the presence of the von-Schweidler-law process as a 
universal feature (cf. Sec. HVDj) . such a polymer-specific feature shows up after the end of 
the von-Schweidler-law process but before the onset of final diffusion [ijj]. 

In Fig. [9] the data are plotted versus Dt so that the simulated and theoretical curves 
coincide in the diffusive late-a regime. This representation facilitates the comparison in the 
j3 and early-a regimes. (Plotting MSDs versus t/r q * as in Figs. H] and [7] leads to a horizontal 
shift of the theoretical curves to the right by a factor of about 0.3 on the log 10 t axis.) 
Figure [9] demonstrates that the theory in terms of the MCT a-master curve (solid line) 
describes the polymer-specific subdiffusive variation of the a-process, where gu(t) ~ t x with 
x = 0.63. Solving the MCT equations for e MCT = —0.046, whose results are drawn as dashed 
lines in Fig. [9j the description of the simulated gu(t) can be extended to about 7 decades in 
t. Our theory derives the Rouse result, x — 1/2, for — » oo (cf. Appendix IB1 and Ref. [15]), 
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in agreement with simulations of long chains in which topological constraints are eliminated 



341 ] . Thus, we suggest that the somewhat larger x = 0.63 found in our N = 10 model is a 



deviation from the pure Rouse behavior due to finite- iV effects (cf. Appendix IB 4[) . 

The inset of Fig. compares the ratio gi(t)/g^(t) of the MSDs for the end [<7i(£)] and 
central [<7s(£)] monomers in a chain. This ratio is 1 in the ballistic and diffusive regimes, and 
exhibits a maximum for times, where conformational motion dominates the dynamics [i.e., 
for t where gu{t) ~ t ' 63 ]. The ratio is also close to 1 in the (3 regime, indicating that the 
cage effect slows the motion of end and inner monomers in the same way. This is because 
the dynamics in the (3 regime is dominated by the confined dynamics of the monomers 
inside the cage. The inset of Fig. M also reveals semiquantitative agreement between theory 
and simulation for later times where the ratio exhibits a maximum. This indicates that 
the motion of the end monomer and that of the central monomer are well discriminated 
by the theory. Thus, chain-end effects for the dynamics of a tagged chain are properly 
taken into account by our theory, even though this effect was neglected in the (static) direct 
correlation functions [cf. Eq. (fl5l) ]. This is because the matrix structure of Eqs. ()3Tj) and 
( |39l) is preserved for describing the single-chain dynamics. The found agreement of the ratio 
gi{t)/g§{t) also shows that its maximum value, which is somewhat smaller than 2 - the 
result expected from the Rouse theory - reflects deviations due to finite- N effects. 

Concerning gc(t), on the other hand, the theory is not so satisfactory: besides the un- 
derestimated plateau height, a careful examination of Fig. |9] indicates that the theoretical 
gc{t) enters the diffusion regime earlier than the simulated one. We will come back to this 
point in the following. We only notice here that the disagreement in gc(t) does not carry 
over to <7m(^)> because gmif) 3> gc(t) for times before the onset of the final diffusion regime. 

As discussed above, our theory yields a subdiffusive, Rouse-like behavior close to T c . 
Clearly, this polymer-specific feature is also present in the simulation at high T. However, 
as T increases, the cage effect loses its importance, and it is thus a priori not clear to what 
extent the theory can still be applied. To examine this, we analyze in Fig. [10] the MSDs at 
T = 1, which is more than twice T C MD . Here, the theory utilizes S(q) taken directly from 
the simulation at this temperature. Figure [10] indicates that, beyond the short-time regime, 
the agreement between theory and simulation is very good for gu{t)- In particular, we find 
gu(t) ~ ^°' 63 with the same exponent. Thus, though originally developed to describe glassy 
dynamics, our theory can also properly deal with the conformational dynamics in normal 
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liquid states. 

The upper inset of Fig. [TOl exhibits the function guc(t) = 9u(t) — Qc{t)- According to 
Eq. (I49p . this function highlights the contributions from the chain conformational fluctua- 
tions to gyi{t), and depends only on the Rouse modes of nonzero mode indices, p > 1. The 
inset clearly indicates that the subdiffusive behavior ~ t om entirely comes from the Rouse 
modes, and that the exponent 0.63 does not reflect a crossover effect from the pure Rouse 
behavior to the final diffusion (i.e., At 0,5 + 6Dt ^ t om ), but is indeed due to finite- N effects. 

On the other hand, we find at T = 1 again the same disagreement for gc{t) and for the 
ratio T q /r q *: (i) The theoretical gc{t) enters the diffusion regime earlier than the simulated 
one which additionally exhibits a subdiffusive behavior, gc{t) ~ t y with y pa 0.8, known 
as "anomalous CM diffusion" [?], [^j]. (ii) The theoretical a- relaxation times r q of 4>(q,t) 
agree quantitatively with the simulation for q > q*, but not for q pa q c (cf. the lower inset 
of Fig. [10]). Thus, the disagreements observed at T = 0.47 are already present at high T, 
suggesting that they are not directly related to the glass transition. 



VI. SUMMARY AND CONCLUDING REMARKS 



In this paper, we proposed a unified first-principles description of the collective struc- 
tural slowing-down and of the single-chain conformational fluctuations in a melt of unen- 
tangled polymers. The description requires static input which can be taken directly from 
simulations, uses approximations like the equivalent-site approximation that can be tested 



explicitly 



161 ]. and attains semiquantitative agreement with simulation results concerning 



collective as well as single-chain dynamics. Our comparative study of theoretical predic- 
tions and simulation data identifies local structural correlations of monomers as the origin 
for the onset of glassy slow dynamics. It is also shown that the chain connectivity causes 
the polymer-specific long-time anomalies of the a process, which manifest themselves in the 
subdiffusive monomer mean-square displacement. Thus, the widely used picture of polymer 
transport in unentangled melts, the Rouse model - including deviations due to finite N - 
emerges from our first-principles approach (cf. Appendix IB"]) . 

On the other hand, we also found deviations between theoretical and simulation results. 
Though probably not directly related to the glass transition (cf. the end of Sec. IVDp . the 
most noticeable disagreement occurs in the collective density fluctuations on the length scale 
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of a chain, q ~ qc (Fig- E]). This disagreement might explain why our theory underestimates 
T c (T C MCT « 0.277 < T C MD « 0.45). There are discernible slow modes at q ~ gc; not 
accounted for by our theory, which appear to couple to the relaxation at other wave vectors. 
The only way in which our theory can compensate this additional coupling is by making 
the cage effect stronger, i.e., by increasing the first peak of S(q) via a decrease of T C MCT . 
Another possible source for the underestimation of T C MCT could result from the neglect of 
the triple direct correlation functions C3 in our theoretical calculations. It was found from 
an MCT analysis for a model of orthoterphenyl that including C3 considerably increases 



T. 



MCT 



36]. For our polymer model, C3 has been determined from simulations in Ref. 16], 



and in principle it would be rather straightforward to take them into account in our theory. 
Unfortunately, the statistical accuracy of the simulated C3 was not sufficient to allow for a 
meaningful test to investigate to what extent the inclusion of C3 affects the value of T C MCT . 

Inspection of Fig. [10] implies that the disagreement found in the theoretical predictions 
- the one in the collective dynamics at q ~ qc and the other in gc(t) concerning the 
anomalous CM dynamics - might be somehow related, since these are the features for which 
our theory does not work well. Superficially, this conjecture agrees with the physics discussed 



in Ref. 



35J. There, the anomalous CM MSD is connected to the polymer coils interacting 



as spheres of radius of gyration R g , and the dynamics at qc ~ 2ir/Rg reflects the polymer 
packing. This implies that taking into account the spatial correlation of CMs through the 
CM structure factor Sc{q) might improve the theoretical results on <p(q,t) at q ~ qc and 
gc(t). Implementing this idea is rather straightforward (cf. Ref. {37I for a related problem). 
However, no improvement was obtained in our case, certainly because Sc{q) at qc is already 
close to 1 (cf. the inset of Fig. Thus, the static coupling between the CMs in our model 
is very weak. Furthermore, also their dynamic coupling is found to be weak, as evidenced 
by the close agreement of their coherent and incoherent intermediate scattering functions at 
. At present, it is not clear how to improve the theory to account for the deviations 



qc 



observed in the collective dynamics at q « qc and in the CM MSD gcif)- It would be 
interesting to investigate to what extent such features are universal or model dependent. 
Only comparison with other models can elucidate this point. 

There is another interesting related issue concerning the collective dynamics at q ~ qc- 
One observes from comparing Fig. [3]^b) with the lower inset of Fig. [10] that the simulation 
result for the ratio T qc /r q * of the ^-relaxation time at qc to the one at the structure factor 
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peak position q* decreases with decreasing T toward T C MD . A similar feature was observed in 
a simulation result for a model of orthoterphenyl which also exhibits some unusual properties 
at intermediate wave numbers corresponding to qc of the present polymer model y. a 
similar T dependence of the ratio T q /r q * at intermediate q range (~ 0.4?*) was also found 



in the coherent neutron-scattering results for a real polymer system [39] . As discussed in 
some detail in Ref. j^ ]. such T dependence of the ratio of the a-relaxation times is beyond 
the implication of MCT. Thus, further investigations are necessary for a comprehensive 
understanding of the yet theoretically unexplained dynamics at intermediate wave numbers 
which are observable in simulation and experimental data for polymer systems. 



Acknowledgments 

Financial support by the DFG and MENRT (IRTG "Soft Matter"), the IUF, the DAAD 
(Grant No. D/00/07994), the ESF SUPERNET Programme, and Grant-in-Aids for scientific 
research from the Ministry of Education, Culture, Sports, Science and Technology of Japan 
(Grant No. 17740282) is gratefully acknowledged. 



APPENDIX A: DERIVATION OF THE MCT EQUATIONS OF MOTION 

This appendix is devoted to the derivation of the MCT equations of motion for gen- 
eral flexible (in the sense that constituent atoms are bonded by some non-rigid potential) 
molecules. Additional approximations, introduced specifically for handling polymeric sys- 
tems, are discussed in the main text. 



1. Zwanzig-Mori equation of motion 



We start from the derivation of an exact equation of motion for the site-site density 
correlators F q h (q,t) = (p a (q)*e J ' c V&( ( l))/ n based on the Zwanzig-Mori projection-operator 



formalism 



21 



. Here, C denotes the Liouville operator 



n N 



N 



iC 



L ' L d i Q Y a m ' J ' J 



i=l a=l 



m 



i,j=l a,b=l 



<9v, c 



(Al) 
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where r" (v°) denotes the position (velocity) of the site a in the zth molecule. The interaction 
potential U(r) comprises both the intra- and inter-molecular contributions (cf. Sec. HT|) . 

To derive the exact equation for F ab (q,t), let us also introduce the longitudinal current 
density fluctuations j a (q) = Y17=i v tz e * q '** where the wave vector q is chosen along the 
z axis, and v1 z (t) denotes the z component of the velocity. p a (q) and j a (q) satisfy the 
continuity equation 

p a (q) = iCp a (q) = iqj a (q), (A2) 

in which the dot denotes the time derivative. For flexible-molecule systems whose kinetic 
energy reads ^ m(v") 2 /2, the static longitudinal current correlation function is given by 

J ab (q) = -(j a (q)*j'6(q)) = 5 ab v 2 , (A3) 
n 



401 ]. there is no off-diagonal 



due to the equipartition theorem. Unlike for rigid molecules 
element and no wave-number dependence in J a b(q)- 

Let us introduce two row vectors p(q) and j(q), whose components are p a (q) and j a (q)) 
respectively. Combining p(q) and j(q) to form a new row vector C(q) = ( p(q) j(q) ), we 
introduce the following projection operator V which acts on some row vector X(q): 

VX = C(C,C)- 1 (C,X). (A4) 

Here, the inner product of two row vectors, Ai and A 2 , shall be defined as the canonical 
ensemble average, (A l5 A 2 ) = (A\A 2 )/n, in which A\ denotes a column vector adjoint to 
Ai, and the factor 1/n is a matter of convention. The matrix (C, C) can thus be expressed 
in terms of the site-site static correlation functions as 

(CChN ° V (A5) 
V 3(q) J 

and its inverse is trivially given in terms of S _1 (g) and J~ 1 (q). 

With the projection operator V so defined, the standard procedure of the Zwanzig-Mori 
formalism leads to the following equation of motion for F a h(q,t): 

N N t 

F ab (q, t) + J2 Fxb{q, t) + ^ / dti M ax (q, t - f) F xb (q, t') = 0. (A6) 

x=l x=l ^° 

Here the characteristic frequency matrix reads 

= q 2 v 2 S^(q), (A7) 
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and the formal expression for the site-site memory kernel M a t,(q,t) is given by 



M ab (q,t) = —{fa(q)*exp(iQCQt)f b (q)), 



(A8) 



in terms of the fluctuating random force 



/ (q) = i'a(q) - «5U 2 J] Pv(<l) S vM> 



(A9) 



1/ 



which evolves with the generator <2£<2, where Q = 1 — V. So far, no approximation has 
been invoked, and the above equation for F a j,(q,t) is formally exact. 

2. Mode-coupling approximation 

The basic idea behind the mode-coupling theory is that the fluctuation of a given dy- 
namical variable decays, at intermediate and long times, predominantly into pairs of hy- 
drodynamic modes associated with quasi-conserved dynamical variables. It is reasonable to 
expect that the decay of the memory function at intermediate and long times is dominated 
by those mode correlations which have the longest relaxation times. The sluggishness of 
the structural relaxation processes in glass-forming systems suggests that the slow decay of 
the memory function at long times is basically due to couplings to wave-vector-dependent 
pair density modes of the form A Xll (k, p) = p A (k)p^(p). The simplest way to extract such 
slowly-decaying part is to introduce another projection operator V2 which projects any vari- 
able onto the subspace spanned by A\n(k, p). Translational invariance of the system implies 
that the only Ax^k, p), whose inner products with a dynamical variable X(q) are non-zero, 
are for the wave vectors p satisfying p = q — k. From here on, we denote by A\^ those 
A Xll (k, p) in which p = q — k, and we define 



Here the factor 1/2 is to avoid the double counting in the summation over the wave vectors, 
and the inverse is defined via 




(A10) 




(AH) 



It is readily verified that P 2 is idempotent and Hermitian. 
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The first approximation in the mode- coupling approach thus corresponds to replacing 
the time-evolution operator exp(iQ£Qt) by its projection on the subspace spanned by A\^. 
exp(iQ£Qt) « V2exp(iQ£Qt)V2- Under this approximation, the memory function reads 

M Mt) = ^ W«(q)* exp(tQ£Qt)V 2 f b (q)). (A12) 

The second approximation is to factorize averages of products, evolving in time with the gen- 
erator Q£Q, into products of averages formed with variables evolving with £ (factorization 
approximation) : 

^(p A (k)>>)V 2£2 Vv(k)p M '(p)> « F xx ,(k,t)F^,(p,t). (A13) 

Specializing this approximation to t = 0, it follows from Eq. (lAlip that the denominator in 
Eq. (IA10I) is given by 

(A x „ Ay,,)- 1 = \ S$(k) Sft(p). (A14) 
Let us obtain the explicit expression for the projected random force, 

V 2 fM) = P 2 j Q (q) - nv 2 ^2(P 2 p v (q))S^(q). (A15) 

V 

To this end, we need to evaluate triple correlations (A Xfl , Ja(o)) an d (A^, Pv(s\))- The former 
can be expressed as 

(pA(k)p„(p),j a (q)) =^ z -(^(k)p;(p)j a (q)}+^i(p*(k)i*(p) Ja (q)), (A16) 

where we have used the relation (AB) = —(AB) and the continuity equation (jA2j) . fc z 
(p 2 ) denotes the z component of the vector k (p). Since (v£ z Vj Z ) = 6ij6 a bV 2 , there hold 
{fxO^P^i&jai^/n = 8 a xv 2 S Xl ,(p) and (p A (k)j*(p)j a (q))/n = 5 a „v 2 S Xfi (k), leading to 

(p A (k)p M (p), ja(q)) = ik z 5 aX v 2 S\,(p) + ip z S a ,v 2 S Xfl (k). (A17) 

The other triple correlation can be expressed in terms of the three-site static structure factor: 

(p A (k)/v(p), p„(q)) = -(pA(k)*p At (p)*p I ,(q)) = S A ^(k, p, q). (A18) 



In the present study, the convolution approximation developed in Ref. 40j shall be employed: 



S X/JiV (k, P, q) ~ S Xa {k) Sftaip) S ua (q). (A19) 
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Using these results, we finally obtain 



n 

k 



EE^[ c Aa r W + c Aa (A;)]pA(k)p (p). (A20) 



Here, c^ tra (g) = [5 a & — w db{o)]/p denotes the intra-molecular direct correlation function 41], 
whereas c a b(q) is the inter-molecular one defined in Eq. (fT4l) . Thus, the projected random 
force naturally comprises the intra-molecular as well as inter-molecular contributions. 

The MCT expression for M ab (q,t) can be obtained by substituting Eq. flA2Dj) into 
Eq. (1A12I) and then invoking the factorization approximation (1A13j) . but let us make here 
a comment on the intra-molecular contribution. Within the same approach outlined above, 
one can derive the MCT expression for the memory kernel for the site-site transversal current 
density correlator, whose q — > limit is related to the shear-stress autocorrelation function 
G(t) [2l|]. It is well known for unentangled polymer chains that G(t) exhibits a power-law 
decay ~ t -1 / 2 for long times. According to the Rouse theory, this polymer-specific decay 
in G(t) is accounted for by the intrachain (or Rouse-mode) contributions However, we 
found that our intra-molecular contributions given in terms of c intra do not lead to this Rouse 
model result for G(t). This implies that a completely different kind of approach is necessary 
for a proper treatment of the intra-molecular contributions in the coherent moduli. Indeed, 
we found a reasonable approach starting from a different projection operator for the intra- 
molecular forces, which reproduces the Rouse model result for G(t). This issue, however, 
shall not be investigated further here, and will be studied in a forthcoming article. Let us 
only mention that (i) even with the inclusion of such intra-molecular contributions to the 
random force or to the memory kernel, we confirmed that all the theoretical results presented 
in the main text are not much affected, and (ii) neglecting intra-molecular contributions to 
the fluctuating force does not mean that intra-molecular couplings are completely discarded 
in our theory, since the intrachain static correlations are properly taken into account via 

Wab{q)- 

The following expression for V^faio) shall therefore be employed in the present work: 

V 2 f a (q) = E E@ • PA(k)p.(p), (A21) 

k A 

in which we have expressed k z as q • k with q = q/q. With the use of the factorization 
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approximation (1A13j) . this leads to 

2 N r 

M ab (q,t) = J dk{(q-k) 2 c Xa (k)c IMb (k)F^(k,t)F ab (p,t) + 



+ (q- k)(q • p)c Xa (k)c^ b (p)F Xb (k,t)F ail ( P ,t)}. (A22) 



3. MCT equations for tagged molecule's correlator 

The MCT equations of motion for a tagged molecule (labeled s) can be derived in a similar 
manner, and only the outline of the derivation and resulting equations shall be presented 
in the following. The Zwanzig-Mori equation for the tagged molecule's density correlator 
F^ b (q,t) = (Pa(q)e l£ 'pfe(q)) is obtained from the projection operator V s onto p s a (q) and 
jKq) = t '^ e * q r ") an d is given by 

N N t 

K b (Q,t)+J2 n ^ F ^) + J2 / dt'M s ax ( q ,t-t')F s xb (q,t')=0. (A23) 

x=l x=l ^° 

Here the characteristic frequency matrix reads 

= q 2 v 2 w- ab \q\ (A24) 

and the formal expression for the memory kernel is given by 

M s ab (q,t) = l(^(q)*e X p(iQ s £Q s t)/ 6 s (q)}, (A25) 
v 

in terms of the fluctuating random force evolving with Q S CQ S (Q s = 1 — V s ) 

/ a s (q) = £(q) - nv 2 p 8 M^M- (A26) 

V 

The memory kernel under the mode-coupling approximation reads 

M s ab (q,t) = \(Vlf:(qyexp(zQS£QH)Vyi(q)). (A27) 

v 

Here the operator V\ projects any variable onto the subspace spanned by pair density modes 
A\^(k, p) = p^(k)p A1 (p) formed by the tagged molecule's density fluctuations and collective 
ones. Adopting the convolution approximation relevant here |40[] 

(p s A (k)*p M (p)*pt(q)) «^^WpV(p)^(?). ( A28 ) 
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one obtains for the projected random force 



k X,/i,u 

2 E E {M^W [I " pc intra (p)] a - 1 c„(p)pl(k)p,(p),(A29) 



pv 
-i — 

n 

k X,fj,,u 



where in the final equality the intramolecular structure factor is expressed in terms of c intra . 
Again, intra-molecular contributions given in terms of c intra shall be neglected here, which 
leads to 

^ 2 7 a s (q) = E EG • P) c <^) AMpM- (A30) 

k A 

All the theoretical results presented in the main text are not much affected by this neglect 
since w ab (q) for microscopic wave vectors, in particular, near the peak position q* of the 
static structure factor S(q), are close to diagonal [if]], in which case there holds c^ tra (g) ~ 0. 
In addition, our derivation of the Rouse model presented in Appendix [B] is not altered since 
the relevant memory kernel there, rh\=o(t), is formed by the summation Yla=i ^faitf ~ > 0) 
of the projected random forces in the small wave vector limit, for which both Eqs. f)A29j) 
and ( 1A30I) yield the identical expression. Under the factorization approximation 

^<P:(k)VA(p)V ss ^ s Vs(k)p,(p)> » F s ab (k,t) F x ,(p,t), (A31) 

substituting Eq. fTA"30]) into Eq. flA"27|) finally yields 

2 N r 

M s ab (q,t) = E J dk^.p) 2 c aX (p)cUp)F s ab (k,t)F x ,(p,t). (A32) 

APPENDIX B: DERIVATION OF THE ROUSE MODEL 

In this appendix, we show that our microscopic formulation for the polymer dynamics 
based on MCT reduces to the Rouse model in the asymptotic limit of large degrees of 
polymerization N. Implications of our theory in this limit on the Rouse-mode correlators in 
the f3 and a relaxation regimes and possible finite- A^ corrections are also discussed. 

1. MCT equations for mean-square displacements 

We start from deriving the MCT equations for monomer MSDs 

ArM = 5rl b {t) - (Sr^(O) with Sr^t) = <[r a (t) - r 6 (0) ] 2 >. (Bl) 
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Notice that the CM MSD gc(t) and the monomer-averaged one <?m(£) introduced in Sec. lIVCl 
can be expressed in terms of Ar 2 b (t) as 

1 N 

9o(t) = W2 E (B2) 

a,6=l 

1 - 1 

9mW = ^E Ar -(*) = ^Tr[Ar 2 (t)]. (B3) 

a=l 

Since F a s fc (g -»• 0,i) = 1 - q 2 5r 2 ab {t)/Q + 0(g 4 ) [cf. Eq. (I36) ]. the Zwanzig-Mori equation for 
Ar 2 b (t) can be derived from the small-g behavior of Eq. flBTj) 

w n t 

-d'Arl(t)+q 2 J2^(Q^0)Ar 2 xb (t)+Y^ / ^m M (t - t') d t ,Ar 2 xb (t') = 6 5 ab , (B4) 

with the memory kernel m a b(t) = \im g ^ q 2 J2 x w ax(l) m xb(t) [cf Eqs. (|39|) and (140]) ]: 

(*) = ^J dk k 4 S(k)c(k) 2 F s ab (k, t) (f>{k, t). (B5) 

In Eq. (1B4j) . the term q 2 w~ b (q — > 0) has to be kept since w a b(q = 0) = 1 is singular and its 
inverse does not exist. For gc{t), a somewhat simplified equation can be derived by taking 
( W 2 ) E a , b of EQ- (USD and noticing lining 2 = 0: 

-2 d hdt) + — 2 dt'm ax (t-t')d t ,Ar 2 xb (t>) = -. (B6) 

a,x,b=l " 

2. Derivation of the Rouse model as asymptotic solution 

In general no simplification of the complicated couplings in the motions of all monomers 
is possible, as an exact diagonalization of the integro-differential equations for the matrix 
Ar 2 b (t) is required. Only for long times and large degrees of polymerization N, an asymp- 
totic solution is possible and is presented in the following. It rests upon the property of 
the memory functions m a b(t) in Eq. flB5j) that they are 'cut-off' by the collective density 
fluctuations. The slowest collective correlator <p(q, t) is connected with the average monomer 
separation and lies at the position q* of the peak of S(q) (cf. Sec. IV Aj) . Thus, at long times, 
the tagged polymer's density correlator F^ b (q, t) at the (asymptotically iV-independent) wave 
vector q* dominates the memory functions. It is bounded by the intrachain structure factor 
at that wave vector, i.e., w a b(q*) > F* b (q*,t) [in the sense that w a b(q*) — F* b (q*,t) is positive 
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definite]. As long as w a b{q*) on this length scale contains no anomalous correlations extend- 
ing over large monomer separations, i.e., w a b(q*) — > for | a — b | — > oo, the same property 
holds for the memory function as well: m ab (t) — > for | a — b \ — > oo. This property, and 
that m ab (t) decays to zero for times longer than r 9 » , the structural relaxation time, are the 
central ingredients to the derivation of the Rouse model within our approach. Note that 
the above reasoning also holds if the role of the microscopic wave vector q* is replaced by 
some other wave vector q\ as long as q^ is characteristic of local motion and asymptotically 
iV-independent. Thus, the following reasoning also applies to high temperature, where, as 
discussed in Sec. IV D\ the slowest collective mode in the simulation result is found to lie at 



q ~ qc 



42]. 



For long times, t — > oo, Eq. (1B6|) is solved by a uniform increase of all MSDs following 
the CM motion: 

g c (t)^QDt and Sr 2 ab (t) - g c (t) + 0(t x ), (B7) 

where the Markovian limit in the memory function requires t ^> r 9 ». Below, we will deter- 
mine the leading correction that exhibits a power-law behavior with the exponent x = 1/2. 
As explained above, the site dependence of m ab {t) for long times is dominated by the tagged 
polymer's density fluctuations at microscopic wave vectors. Therefore, the summation over 
site indices and consecutive fc-integration will asymptotically become iV-independent, and 
the diffusion constant D scales like 

D = ^r (tf-oo), (B8) 
with the (asymptotically ^-independent) friction coefficient ( determined by 



71 Jo J L 6=1 



(B9) 



15l | follows the asymptotic law (IB8j) 



We confirmed that the Gaussian chains studied in Ref. 
for large N. 

In considering internal- mode contributions to the monomer MSD in the limit of iV — > oo, 
chain-end effects can be neglected, and the structure of the various matrices in Eq. flB4j) 
simplifies. We can assume that they depend only on the difference of indices, s = a — b 
(characteristic of 'Toeplitz' matrices |43j), and we define, e.g., ArJ s=a _ b ^(t) = Ar^(t). This 
assumption neglects monomer correlations caused by chain ends and does not hold, e.g., 
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for isolated self-avoiding-walk polymers in good solvents whose end regions are slightly less 
swollen than middle portions. Similar correlations have recently been discovered also in long 
chain polymer melts, but the amplitude is much weaker than in dilute solution {44]. Here, 
as a first step, we assume that such non-trivial correlations to be absent. 

For the monomer-averaged MSD, it suffices to obtain the distribution of the internal 
modes since gu(t) m Eq. (1B3|) is expressed as the trace of the matrix Ar 2 b (t), where only 
the eigenvalues of the internal modes enter. The distribution in the N — > 00 limit can 
be found by assuming periodic boundary conditions and performing a Fourier transform. 
Transformed quantities like 

00 

= £ eiXsAr U^ ( B1 °) 

s=— OO 

- — -2 

shall be marked by a hat. Notice that there hold gc[t) = (l/iV)Ar A=0 (t) and <7mC0 = 
— 2 

(1/iV) 2jj> Ar x (t), so that the internal-mode contribution to the monomer-averaged MSD 
is given by gyiit) — gcif) — Sa^o ^"aW- m ^ ne asymptotic N — > 00 limit, the 

monomer- averaged MSD follows from the density of states of internal modes via [43 ] 

J^ Ar A(*)" ( B11 ) 

Here and in the following, A 7^ shall be assumed unless stated otherwise. 

2 

The equation of motion for Ar x (t) is obtained from Eq. (1B4I) via Fourier transforma- 
tion, recognizing that matrix products, owing to the assumption of the dependence on the 
index-difference only, become convolution and turn into simple products after Fourier trans- 
formation: 

d 2 t Ar\{t) + r x &l(t) + / dt'm x (t - t>) d t ,Ar\{t') = 6. (B12) 

Jo 



—dfAr 

tr 

Here we have introduced 



f^S^'lS^* (B13) 



s=— oo 



and the transformed memory kernel is given by 

Pi 



m x (t) = dk k 4 S(ky(k)F s x (k, t)(j)(k, t). (B14) 

67T J 

Equations (1B11I) and (1B12I) yield the iV-independent growth of the monomer-averaged 
MSD resulting from the internal modes, whose spectrum shall be determined to lowest 
order in the mode parameter A. The Gaussian approximation shall be assumed for the large 
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separation behavior in Sr^ s=a _ b ^(0) 44j, so that the small A properties of r A can be found 
from 

i 00 2 

f 7 1 « -- Y e as \s\a 2 s ^%. (B15) 
6 3A v ; 

Here <r s denotes the statistical segment length. At the same time, the memory kernel is 

Taylor-expanded for small A, and then a Markovian approximation (cf. the next subsection) 

is performed 



rh\(t) 







8(t) = ((/k B T)8(t), (B16) 



dt m A=0 (t) 

to derive the long time behavior. It is justified, as discussed above, because the memory 
kernel is dominated by microscopic wave vectors, where the summation over s = a — b 
converges rapidly and an expansion in A is possible. In Eq. (1B16j) only the lowest order in 
A is retained, which is given by the friction coefficient from Eq. (1B9I) of the CM motion. 
Neglecting the inertia term for long times, we find from Eq. (1B12j) for the small-A modes 

C , _ 2 3A 2 - -2 

-±-d t Ar x (t) + — Ar A (t) = 6, (B17) 

Kb-L cr s 

with the initial value Ar A (0) = 0, or equivalently 

2 



1 B 

whose solution reads 



Aj^r A (*) + ± r*r A (t) = 0, (B18) 
kuT erf 



5r\{t) = e - 3k » Tx2t ^ <fr A (0), (B19) 

^ 2 

with 5r A (0) = — 2<jg/A 2 . From this, the following monomer-averaged MSD follows for long 
times as familiar in the Rouse model 21: 



9M(t) _ 9c(i) . [J. i (1 _ e -3^) . M ^ ^ ,B 2 0) 

This concludes the derivation of the Rouse model as the asymptotic large chain-length limit 
of the MCT equations for a polymer chain dissolved in a melt of identical polymers. (The 
use of the Markovian approximation will be justified in the next subsection.) We find the 
expected scaling of the diffusion coefficient with molecular weight in Eq. (1B8I) . the (low- 
lying) spectrum of eigenvalues in Eq. (lB19p . and the resulting anomaly in the monomer 
MSD, Eq. flB20p . The occurring parameters can be measured from global chain properties: 
the friction coefficient £ from the averaged friction kernel, Eq. (1B9I) . and the segment length 
a s from the Gaussian behavior of the equilibrium segment correlations at large separation, 
Eq. flBTol) . 
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3. Implications on the Rouse-mode correlators 

^ 2 

We notice that Sr x (t) is essentially the (diagonal) Rouse-mode correlators C pp {t) intro- 
duced in Sec. II V CI with the correspondence A ~ p/N since there holds 6^ p>1 [C pp (0) — 
Cp P (t)} = J2x^ol^ r xi^) ~$ r \(,Q)} from the comparison of Eqs. (14"§|) and (1B11[) . The difference 

^ 2 

comes from the boundary condition adopted in defining C pp (t) and Sr x (t). So, the small-p 
properties of the normalized Rouse-mode correlators c p (t) = C pp (t)/C pp (0) can be deduced 
from the small-A behavior of 

c x (t) = 5rl(t) / 5rl(0). (B21) 
The equation of motion for C\(t) can be derived from Eq. (1B12I) : 
1 ~ f 

-<9 2 c A (t) + T A c A (t) + / dt'm x {t - f) d t ,c x (f) = 0. (B22) 
v Jo 

The Laplace transform of this equation reads 

- —z[l + zc x {z)} + f x c x {z) - m x {z) [1 + zc x {z)} = 0, (B23) 
v 

where the convention f(z) = i f^°dte lzt f(t) with Imz > is adopted. 

Let us consider liquid states, for which there is no nonergodicity pole in the Laplace 
transform of correlators. (Nonergodicity parameters shall be discussed below.) Then, we 
have from the z — > limit of Eq. (1B23I) : T\C\(z — > 0) = rh\(z — ■> 0). Retaining only the 
leading-order contribution in A for T x and m x , one obtains for small A 

2 

cx(z - 0) = ^^x=o(z 0). (B24) 
Since c\(z —> 0) and m x= o(z — > 0) are proportional to their relaxation times, this implies 

^ 2 

that the relaxation time of c{t) [and hence of Sr x (t)} is larger by a factor of 1/A 2 ~ {N/p) 2 
than that of rh x= o(t). This justifies the use of the Markovian approximation, which has been 
adopted in Eq. (1B16j) . Since, as discussed above, the relaxation time of m x=0 (t) is dictated 
by that of the coherent dynamics at microscopic wave vectors, Eq. (IB24I) also implies the 
separation of the time scale for c(t) from that for the density fluctuations at the monomer 
length scale. 

We next turn our attention to nonergodicity parameters f\ = c\(t — > oo) and m x = 
m x (t — > oo) in glass states. Since there holds lim t ^ 00 /(t) = — lim z ^ zf(z), one obtains 
from the z — > limit of Eq. (IB23I) : f x = m x /[m x + T x \. Because of Eq. flB15j) . we have 
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fx = 1 in leading order for small A. This holds also at the MCT critical point, so that the 
critical nonergodicity parameter for c A (i) is = 1. Since c A (i) < 1, this also implies that 
the critical amplitude for c x (t) is zero, h x = 0. Therefore, our theory predicts that, for small 
mode indices or in the N —>■ oo limit, the critical nonergodicity parameter of the Rouse-mode 
correlators is unity, and they do not exhibit the MCT j3 dynamics (cf. Sec. IIVDI) . 

Let us then consider the dynamics of c\(i) in the a regime, which can be described by its 
a master curve [cf. Eq. (1601) ] . We consider the dynamics on the a- relaxation time scale t' a 
[cf. Eq. (1571) ] and write, e.g., c x (t) = c\(t) with t = t/t' a and c x (z) = t' a c x (z) with z = zt' a , 
with a yet unspecified function c A . It then follows from Eq. (]B23[) : 

- \ [1 + zc x {z)\ + f x c x (z) - rh x {z) [1 + zc x {z)\ = 0. (B25) 

Now, the a-scaling limit shall be performed: t' a — > oo for T — ► T c +, but with t and z 



fixed [11]. We thus obtain T x c x (z) — m x (z) [1 + zc\(z)\ = 0. Since c x {i — > 0) = f^ 11] and 
fl = 1 for small A as derived above, the inverse Laplace transform of this equation yields 

T x c x (t) + f £ m x (t - ?)^c A (?) = 0. (B26) 
J o 

Again, only the leading-order contribution in A for T x and shall be retained. With the 
same reasoning as presented concerning Eq. flB24j) . one obtains the time-scale separation of 



the dynamics of c x (t) from that of fh x=0 (t). Therefore, the Markovian approximation for 
Eq. (IB26P is justified, leading for small A to the exponential decay of the a master curve 

c A (*) = exp[-f/f A ], (B27) 

with the relaxation time f x = (a^ / 3/cbTA 2 (£ = limr- > r c + C/O whose dependence on the 
mode index reads 1/A 2 ~ (N/p) 2 . Thus, our theory in the N oo limit yields the Rouse- 
model result for the a-master curves for the Rouse-mode correlators. 



4. Finite-N corrections 

So far, we have derived the asymptotic solution of our MCT equations in the limit of large 
degrees of polymerization N, by retaining only the leading contribution in the expansion 
in the mode parameter. To find finite-iV corrections, one has to go beyond the leading 
order, but it is difficult to explicitly work this out. On the other hand, it is obvious that 
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finite- N corrections lead to deviations from the asymptotic (Rouse model) results: when the 
normalized Rouse-mode correlator c p (t) in the a regime is fitted via a Kohlrausch function 
A p exp[—(t/T p ) /3p ], finite-iV corrections lead to A p < 1, (3 P < 1, and deviations from r p oc 
(N/p) 2 . Furthermore, the monomer-averaged MSD does not exhibit the square-root-time 
dependence any longer. 
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FIG. 1: (Color online) Collective static structure factor S(q) of the melt as a function of the 
modulus of the wave vector q for temperatures T = 0.47 (solid line), 0.70 (dashed line), and 1 
(dotted line). S(q) exhibits a maximum around q* = 6.9 whose position is indicated by an arrow. 
The inset shows Sc(q), the static structure factor of the chain's center of mass, for T = 0.47, 
0.70, and 1. There is practically no temperature dependence in Sc(q), and three curves cannot 
be distinguished from each other. Sc(q) exhibits a weak maximum at qc = 3.4 whose position is 
indicated by an arrow. 
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FIG. 2: (Color online) Comparison of the static structure factor S(q) (circles) with the site- 
dependent static structure factors S a (q) for a = 1 (dashed line), a = 2 (solid line), and a = 5 
(dotted line). The inset compares S(q) (circles) with 1/S~ 1 (q) for a = 1 (dashed line), a = 2 (solid 
line), and a = 5 (dotted line). (The dotted lines for a = 5 in the main panel and in the inset are 
not clearly visible since they almost agree with the solid lines for a = 2.) S a (q) and S~ 1 (q) are 
defined by the first equality of Eqs. (fl9|) and (|20l) . respectively. S(q), S a (q), and 5 , ^ 1 (g) are taken 
from the simulation at T = 0.47. 
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FIG. 3: (Color online) (a) Glass-form factors f c (q) of the coherent density correlators <ft(q, t) versus 
q. The circles represent the result from the simulation at T = 0.47, and the solid line that from 
MCT. The dashed line denotes the extrapolated S{q) (multiplied by 0.1) at T C MCT ps 0.277. The 
arrows indicate the peak positions q* and qc of S(q) and Sc(q) (c/. Fig. [Q. The inset depicts the 
extrapolated S(q) at T c MCT (dashed line), and the simulated S(q) at T = 0.47 (solid line), 0.48 
(dotted line), and 1 (long-dashed line) around the peak q* . (b) Rescaled a-relaxation times r q /T q * 
(main panel) and the stretching exponent (3 q (inset) of 4>(q, t) versus q. The circles represent the 
result from the simulation at T = 0.47, and the solid line that from MCT. 
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FIG. 4: (Color online) (a) (ft(q,t) as a function of t/r q * for q = 4.0 (left scale), 6.9 (right scale), 
and 12.8 (right scale). t 9 * is the a-relaxation time at q*. The circles refer to the simulation results 
at T = 0.47, the solid lines to the MCT a-master curves, and the dashed lines to the MCT curves 
at the distance parameter e MCT = —0.046. (b) 4>(q,t) as a function of t/r q * for q = 6.9 and 12.8. 
The circles and the dashed lines are the same as in (a), but here the dotted lines denoting the 
MCT curves at the distance parameter e MCT = —0.022 are included as well. 
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FIG. 5: (Color online) Glass-form factors f sc (q) of the correlators (j) s (q,t) (a) and f pc (q) of 
the correlators (f) p (q,t) (b) as a function of the wave number q. The circles represent the result 
from the simulation at T = 0.47, and the solid line that from MCT. The dash-dotted line in (a) 



denotes /^(q) based on the Gaussian approximation (|64p with the value 



M 



0.098 taken from 



the theoretical calculation. The dashed line in (b) shows the simulated w(q) (multiplied by 0.1) at 
T = 0.47. 
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FIG. 6: (Color online) Single-chain density correlators <j) s (q,t) (a) and cp p (q,t) (b) as a function 
of t/r q * for q = 4.0, 6.9, and 12.8. r q * is the a-relaxation time of the coherent density correlator 
<f>(q,t) at q = q*. The circles refer to the simulation results at T = 0.47, the solid lines to the MCT 
a-master curves, and the dashed lines to the MCT curves at the distance parameter g MCT = —0.046. 
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FIG. 7: (Color online) (a) Normalized Rouse-mode correlators c p (t) = C pp (t) / C pp (0) as a function 
of t/r q * for p = 1, 2, 3, 5, and 9 (from right to left). r q * is the a-relaxation time of the coherent 
density correlator (j)(q,t) at q = q*. The circles refer to the simulation results at T = 0.47, and 
the solid lines to the MCT a-master curves, (b) Enlargement of the (3 region in (a); the results 
for p = 9 are omitted. Here, dashed lines represent the MCT curves at the distance parameter 

e MCT = _ a04 g_ 
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FIG. 8: (Color online) The plateau heights f£ (a), the ratio T p /r q * of the a-relaxation times (b), 
and the stretching exponent j3 p (c) of the Rouse-mode correlators c p (t) as a function of the mode 
index p. The circles represent the result from the simulation at T = 0.47, and the solid line that 
from MCT. The dotted line in each panel refers to pure Rouse behavior predicted by our theory 
in the asymptotic limit of large N (cf. Appendix [B]) : f p = 1, r p ex [sin(p7r/2A r )]~ 2 , and (3 p = 1. 
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FIG. 9: (Color online) Double logarithmic presentation of the MSDs <?m(£) (labeled M, left scale) 
and gc(t) (labeled C, right scale) as a function of Dt. The inset exhibits the ratio gi(t)/g^(t) (end- 
over middle-monomer MSD). The circles refer to the simulation results at T = 0.47, the solid lines 
to the MCT ce-master curves, and the dashed lines to the MCT curves at the distance parameter 
£ mct _ _o.046. The dash-dotted lines indicate the diffusion law, 6Dt, while the dotted line shows 
the power law, ~ t 0,63 . 
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FIG. 10: (Color online) Double logarithmic presentation of gyi(t) (labeled M) and gc(t) (labeled 
C) versus Dt at T = 1. The upper inset exhibits gMc(t) = gM(t) — gc(t), whereas the lower inset 
shows the q dependence the ratio T q /r q * of the a-relaxation times of the coherent density correlators 
(p(q,t) at T = 1. The circles represent the result from the simulation, and the solid line that from 
MCT. The dotted line in the main panel and the upper inset denotes the power law ~ i - 63 . 
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